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We report here an assessment of carrier multiplication (CM) yields in PbSe and PbS nanocrystals 
(NCs) by a quantitative analysis of biexciton and exciton dynamics in transient photoluminescence 
decays. Interest in CM, the generation of more than one electron and hole in a semiconductor after 
absorption of one photon, has renewed in recent years because of reports suggesting greatly increased 
efficiencies in nanocrystalline materials compared to the bulk form, in which CM was otherwise too 
weak to be of consequence in photovoltaic energy conversion devices. In our PbSe and PbS NC 
samples, however, we estimate using transient photoluminescence that at most 0.25 additional e- 
h pairs are generated per photon even at energies huj > 5E g , instead of the much higher values 
reported in the literature. We argue by comparing NC CM estimates and reported bulk values 
on an absolute energy basis, which we justify as appropriate on physical grounds, that the data 
reported thus far are inconclusive with respect to the importance of nanoscale-specific phenomena 
in the CM process. 



I. INTRODUCTION 

The process of carrier multiplication (CM) consists of 
the generation of more than one electron and hole after 
absorption of a single photon in a semiconductor, its ef- 
fectiveness determined by a rich interplay of the interac- 
tions between charge carriers, phonons and lightP^From 
a practical perspective, though, its chief potential as an 
enabler of more efficient solar spectrum harvesting in en- 
ergy conversion devices has been limited by the very weak 
CM of bulk materialsP The topic of CM has however 
reemerged in recent years due to reports of very strong 
enhancements of the CM process in nanocrystalline semi- 
conductors. 

Enhanced CM was first reported for PbSe and PbS 
nanocrystals (NCs) by Schaller et alP and Ellingson ct 
alP using the transient absorption (TA) technique. Work 
on this material system has been extended, with one re- 
port inferring the creation of up to 7 e-h by a single high 
energy photon based on pump-probe dataP and a study 
suggesting that the enhancement occurs not only for NCs 
in solution but also in close-packed films relevant for po- 
tential device applications. 6 Other material systems have 
also been explored, with work ini tially sho wing evidence 
for strong CM as well in CdSepElnAsP^and Si NCsP 

Since then, there have been several reports observ- 
ing little or no CM. Using a transient photolumines- 
cence experiment, we found no evidence for CM in 
CdSe NCs at energies well above previously reported 
thresholds More recently, Pijpers et al. have re- 
ported difficulty in reproducing their observation of CM 
in InAiP, and a new study has reported no observable 
CM in InAs/CdSe/ZnSe (core/shell/shell) NCsPln ad- 
dition, there remain several unresolved questions pertain- 
ing to CM in lead chalcogenide NCs. For instance, there 
are significant qualitative and quantitative differences be- 
tween the Schaller et al. reports^^ of strong CM follow- 
ing a universal trend with ftio/E g and the Ellingson et 



al. results^ which appear to show a distinct particle size 
dependence and smaller yields, in some cases by factors 
of 2-3. Second, considerable theoretical debate about 
CM in NCs remains, mostly due to a lack of information 
about intraband relaxat ion processe s deep in the exciton 
and biexciton manifolds! 16 * 17 * 18 * 19 * 20 ! Recognizing this de- 
ficiency, Allan and Delerue have allowed for a wide range 
of intraband relaxation rates in their flexible theoretical 
framework, but still find the largest CM yields reported 
by Schaller et al. difficult to accommodate.^ Overall, 
these outstanding issues suggest the need for continuing 
the assessment of CM in lead chalcogenide NC samples. 

In this work we study carrier multiplication in PbSe 
and PbS nanocrystals using transient photoluminescence 
(tPL), a technique that more specifically informs on the 
e-h pair population within NCs t han the pump-probe 
methods commonly employed! 12 * 21 * We first characterize 
the exciton and multiexciton PL signatures in these ma- 
terials using low photon-energy excitation. We find that 
PbSe and PbS NCs, when adequately surface passivated, 
have flat exciton population dynamics over a Ins window. 
At higher excitation power, strong features appear with 
fast 50-200ps decay lifetimes attributed to biexcitons. Af- 
ter these calibration steps, we measured tPL decays to 
look for evidence of CM using excitation at 3.1eV, well 
above previously reported CM energy thresholds for the 
NC materials in this studyP^ Although we distinctly ob- 
serve a signal consistent with CM for all of our PbS and 
PbSe NC samples, the CM yields we estimate, defined 
as the average number of additional e-h pairs generated 
per absorbed photonP^ reach only y cm w 25% even when 
huj > 5E g . These values are significantly lower than those 
of previous reportsP^ 

In the final section we explore the issue of comparing 
CM yields between NCs of different sizes and with the 
bulk. We show that if nanoscale-specific physics, such as 
potentially slowed intraband relaxation, are not a priori 
assumed, one would expect CM yields to depend only on 
the incident photon energy, regardless of the particle's 
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FIG. 1: Absorption spectrum of a typical PbSe NC sample 
used in this study. The NC bandgap E g — 0.84eV is deter- 
mined as the peak of the first absorption feature. 



size. This suggests that CM yields be compared on an 
absolute photon energy basis. Revisiting the literature in 
this framework shows that the reports on PbS and PbSe 
NCs to date do not uniformly suggest very large enhance- 
ments of the underlying CM physics when compared to 
what has been reported^ for bulk PbS. 



II. EXPERIMENTAL 

PbSe and PbS NCs were prepared by high temper- 
ature pyrolysis of Pb andSe/S precursors in an oleic 
acid/octadecene mixture! 23 * 24 ! The growth solutions were 
purified by a single precipitation, redispersed in hexane, 
and transferred to 1mm path length quartz cuvettes in 
a nitrogen glovebox. The resulting samples, with opti- 
cal densities of ~ 1-3 at 1.55eV, were sealed and taken 
out into air for subsequent measurements. As will be 
described below, some samples of larger particles (first 
absorption feature < 0.8eV) were treated with Cd 2+ by 
adding a few drops of cadmium oleate to the hexane NC 
dispersions at room temperature. All samples were mag- 
netically stirred during acquisitions, and PL decays un- 
der weak 1.55eV excitation were periodically monitored 
to check for any degradation. A typical sample's absorp- 
tion spectrum is shown in Fig. [I] 

Transient photoluminescence decays of the samples 
were collected using a fluorescence upconversion appa- 
ratus based on an amplified Ti:sapphire system operat- 
ing at 250kHz. A portion of the pulse train was passed 
through a BBO crystal to generate excitation sources 
at 3.1 eV and 1.55 eV which were separated with two 
dichroic mirrors and focused on the sample to spot sizes 
of roughly ~ 100/im and ~ 50/im diameter respectively. 
Emission was collected in a front-face geometry using off- 
axis parabolic mirrors and focused onto another BBO 



crystal. Following a variable delay, the rest of the 1.55eV 
pulse train was overlapped with the collected emission 
and the resulting sum frequency generation was sepa- 
rated spatially and spectrally using interference filters 
and a monochromator. The signal was detected with a 
cooled PMT and amplified using a lock-in amplifier. For 
these experiments, the pulse width was maintained rel- 
atively long by tweaking the amplifier compressor away 
from its optimal short-pulse configuration to avoid excess 
noise from continuum generation in the mixing crystal. 
We have nevertheless maintained a time resolution bet- 
ter than ~ 15 ps as measured from the rise time of the 
tPL signal. Because the peaks of the exciton and mul- 
tiexciton PL were not found appreciably different within 
our spectral resolution, all decays for a given sample were 
acquired at a fixed wavelength. 



III. RESULTS AND DISCUSSION 

A. Exciton decay dynamics 

We began by characterizing the PL dynamics of single 
excitons (X) in PbS and PbSe NCs using weak 1.55eV 
excitation. In general, samples of small and moderate 
sized NCs had flat PL dynamics over the full temporal 
range of our instrument (see Fig. [2J. In contrast, as- 
prepared CdSe core particles almost invariably show sig- 
nificant sub-na nosec ond dynamics, attributable to trap- 
ping by defects] \l2i2H The flat decays we observe in these 
PbSe and PbS samples suggest goo d sur face passivation 
of NCs prepared by these methods] 23 1 24 1 and are consis- 
tent with the very high luminescence quantum yields re- 
ported in the literature.^ In addition, we also measured 
the PL dynamics over a much longer window for one of 
our PbS samples using an InGaAs amplified photodiode 
and found a nearly single exponential fluorescence de- 
cay with a ~ 660ns lifetime, consistent with previous 
studies .1251 

The PL dynamics of larger as-prepared particles, those 
with Eg <^ 0.8eV, typically showed multiexponential X 
decays with large sub-nanosecond components, suggest- 
ing poor surface passivation. Moreover, these dynamics 
steepened irreversibly upon exposure to 3.1eV radiation. 
In an attempt to remove non-radiative pathways and to 
stabilize the particles, we chose to apply a mild Cd 2+ 
treatment to the NCs. Addition of cadmium oleate to 
hexane dispersions of E g = 0.73eV PbS and E g — 0.60eV 
PbSe NCs resulted in nearly flat single-exponential X de- 
cays and robustness to prolonged 3.1eV irradiation, while 
causing no noticeable changes in the absorption spectra 
and emission wavelength of the samples. Our measure- 
ments also suggest that surface treatment of these sam- 
ples does not have much effect on CM yields. We studied 
one sample of fairly large E g — 0.68eV PbSe NCs that 
did not require Cd 2+ treatment, and found that its es- 
timated CM yields were similar to the other large NC 
samples that were treated with Cd 2+ . We also checked 
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FIG. 2: Transient PL dynamics of a sample of PbSe NCs in 
hexane dispersion (E g — 0.84eV) under weak 1.55eV excita- 
tion. 



the effect of the Cd 2+ treatment by applying it to NC 
samples that already exhibited adequate surface passi- 
vation and found no significant change of the biexciton 
lifetime or estimated CM yield. 

We have chosen to carry out further studies only on 
samples that show flat tPL decays over our experimen- 
tal timescale, whether as-prepared or Cd 2+ -treated, be- 
cause the interpretation of subsequent results is consid- 
erably simplified. A multiexponential X decay entails 
an inhomogcncous distribution of NC surface passivation 
which can then s uppor t a nontrivial distribution of mul- 
tiexciton lifetimes I1212U complicating both the isolation of 
MX features in tPL decays and the quantification of the 
underlying exciton and multiexciton populations. The 
second and more serious problem was that the X de- 
cays of samples with poor surface passivation tended to 
change irreversibly when exposed to 3.1eV for the lengths 
of time necessary to obtain adequately clean data with 
our apparatus. For these reasons we focused only on 
well-passivated samples, ft is conceivable that CM yields 
might depend on the details of the NC surface. If so, the 
results of this work may be difficult to generalize beyond 
the constraints of our particular sample preparation and 
selection methods. 
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FIG. 3: (Color Online), (a) Transient PL dynamics of PbSe 
NCs (Eg = 0.84eV) under increasingly strong 1.55eV exci- 
tation, showing the growth of the BX emission feature (red 
lines). The thin solid black lines are fits to axe 
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clbxc "' with tx > Ins and tbx = 58ps fixed, (b) Fits 
of the X and BX exponential components ax and asx to 
a population profile following poissonian photon absorption 
statistics for an inhomogeneous excitation beam, no denotes 
the maximum average number of photons absorbed. See ap- 
pendix for details. 



B. transient PL of the BX state 

Strong excitation pulses can create biexcitons (BX) 
and further multiexcitons in NCs by sequential photon 
absorption. An excitation power series for our Eg — 
0.84eV PbSe sample is presented in Fig. [3^i, showing 
the growth of a large fast feature, which we attribute 
to the BX, on top of the single X dynamics. These de- 
cays are well described as the sum of a slow X compo- 
nent and a fast BX component, Oflx e ™ + a^e _ */ rx 
with fixed lifetimes tx > Ins and tbx ~ 60ps. Under 



strong excitation, additional faster components appear, 
attributable to emission from higher multiexcitons The 
rapid tbx decay times are due to an Auger-like relax- 
ation mechanisnPS and the rates we measure are consis- 
tent with those measured by pump-probe techniques.^ 

Our method for estimating CM yields,^ described in 
the next section, relics importantly on an accurate cal- 
ibration of the link between observed tPL decays and 
the underlying BX and X populations soon after excita- 
tion. This information can be summarized in the quan- 
tity (clbx l 1 o,x) sat, the ratio of the sizes of the BX and 
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X tPL decay components expected in the hypothetical 
case that all NCs are initially excited to the BX state, 
i.e. when the BX is saturated. In Fig. [3J), we fit the 
observed exponential components, ax and clbx, to pop- 
ulation profiles assuming poissonian photon absorption 
statistics. The power series of X and BX features are 
found consistent with this assumption, and we are able 
to estimate sample-dependent (cibx /ax) sat values in the 
range of 2.5-4. This implies that the radiative rate of the 
biexciton, k r g X , is 3.5 — 5 times greater than k x d . In- 
terestingly, the numbers are similar to those observed for 
CdSe NCs, where we proposed that the enhanced k^ x 
could be due to spin substructure since the lowes t X 
fine structure state in CdSe is known to be darlP^U b u t 
the ground state BX is predicted to be bright™ How- 
ever, it has been suggested that no such spin structure 
is necessary to explain the long X emission lifetimes of 
lead chalcogenide NCs!^ In such a scenario, a simple 
accounting of all the possible electronic configurations of 
band-edge X and BX assuming known selection rules and 
thermal equilibrium gives k^ x — Ak r x d , which is consis- 
tent with our results. A derivation can be found in the 
appendix, along with an explanation of the relationship 
between (asx /ax) sat, k r g X , and k r x d , and a description 
of our population profile modeling. 

C. Carrier Multiplication 

We turn now to studying tPL decays under 3.1eV ex- 
citation. Photons of this energy are well above the CM 
thresholds tha t ha ve been previously reported for PbS 
and PbSe NCs I EE if 

carrier multiplication occurs in our 
samples, it would be reflected in the tPL dynamics as a 
residual BX or higher MX component that persists in the 
limit of very weak excitation, when at most one photon 
is absorbed per NC. Fig. |4^i compares PL dynamics for 
E g = 0.84eV PbSe NCs under 1.55eV and 3.1eV exci- 
tation. As described previously, weak 1.55eV excitation 
results in flat, single exponential dynamics correspond- 
ing to X decay, while at higher power the tPL traces 
exhibit a fast BX component as well. In contrast, even 
at low fluence, excitation at 3.1eV results in decays with 
a fast component closely following BX dynamics. Fig. [4ji 
shows the asx/a-x ratios obtained from a series of mea- 
surements with varying 3.1eV excitation fluence. Our 
extrapolation shows that the BX-like feature persists in 
the zero power limit (P —> 0), and we thus attribute it 
to CM. The CM yield, y cm , for the sample is then given 
byP 



Vcm — 




Because our best estimates of (asx /ax) sat are in the 
range 2.5 — 4, CM yields are smaller by a factor of ~ 3 
than the simple ratio asx/ax would suggest. For this 
E g = 0.84eV sample, y cm w 9% at frw/E g = 3.7. Fig. 
[4}d and [4jl display similar data for a sample of larger 



E g = 0.68eV PbSe NCs. The sample exhibits a bigger 
fast component in the P — > limit, and therefore a larger 
CM yield of w 23% at hu/E g = 4.6. 

We have studied a number of PbS and PbSe NC sam- 
ples in this way and find y cm always in the range of 
10 — 25% even in samples for which fku > 5E g . As is 
summarized in Fig. [5] our CM yield estimates are sig- 
nificantly lower than those previously reported by other 
researchers for their PbS and PbSe NC samples. Schallcr 
et al. report a universal dependence of CM yields on 
hu>/Eg, predictin g CM yields of up to 2 additional e-h 
pairs at Huj = 5-E s P3EHingson et al. do not observe such 
structure in their data, but instead report CM yields that 
appear to depend on particle sizeP The Ellingson et al. 
data fall bellow the universal curve of Schaller at al. by 
a factor of roughly two, while our own best estimates of 
the CM yields are an additional factor of 2-3 smaller. It 
should also be noted that the findings of Schaller et al. 
predict not only BX formation, but also triexciton (TX) 
yields of 0.5 and 1 for our largest samples when excited at 
3.1eV. However, our data fits very well to only a BX and 
an X component. Any appreciable TX would have been 
evidenced in our measured decays since the TX decay dy- 
namics are within our experimental time resolution and 
the TX emission peak is expected to be close to that of 
the X and BX because of the approximate 8-fold degen- 
eracy of the lowest lead chalcogenide NC electron and 
hole states. 

Since our numerical results arc in disagreement with 
the previous reports based on TA techniques, we con- 
sider here possible sources of error in our CM estimates. 
In Fig. [5] we show estimated 95% confidence intervals 
for (aBx /ax)p-^o related to noise in the experimental 
decays. These uncertainties in y cm are all smaller than 
±0.06 and are likely unbiased. The part of our methodol- 
ogy most susceptible to a systematic error is the satura- 
tion ratio of the BX to X tPL components, {asx /ax) sat- 
Any multiplicative error in this quantity translates di- 
rectly into a multiplicative error in the CM yield. In our 
study, we have estimated (a^x I ax) sat by fitting the sizes 
of X and BX decay components under 1.55eV excitation 
to a population profile and then assuming that this satu- 
ration ratio should apply as well to biexcitons created by 
a CM process. Using these (aBx/ a x)sat values of 2.5- 
4, we have determined CM yields in the range of 10-25%. 
For our results to roughly match the magnitudes of CM 
reported by Ellingson et al. we would have had to use 
much smaller values (ogj /ax) sat ~ L with an even fur- 
ther reduction to {asx /ax) sat < 1 required to achieve 
agreement with the Schaller et al. reports. However, us- 
ing such small {a r ^ x / a r x d ) sat would be inconsistent with 
our direct observation of asx > %ax under sufficiently 
strong excitation conditions. We are therefore confident 
in our principal conclusion that the CM yields in the 
PbSe and PbS NC samples we have studied are signifi- 
cantly smaller than those previously reported for the PbX 
material system. 
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FIG. 4: (a) Comparison of PL decays from a sample of E g = 
0.84eV PbSe NCs under 1.55eV and 3.1eV excitation. Even 
as the 3.1eV excitation power reaches the low power limit, the 
decays continue to exhibit a fast component consistent with 
BX dynamics, (b) Same as (a) but for E g = 0.68eV PbSe 
NCs. (c)-(d) Plots of ciBx/ax vs. ax for different weak 3.1eV 
excitation fluences and extrapolation to the ax — > (P — > 0) 
limit for the samples in (a) and (b) respectively. Dividing this 
extrapolated value by the (clbx /ax) sat determined from an 
independent 1.55eV power series gives CM yields y cm = 0.09 
and y cm =0.23 for the two samples at huj — 3.7 E g and huj = 
4.6E g respectively. 
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FIG. 5: (Color online). Summary of CM yields determined 
in this study and comparison to literature reports on PbSe 
and PbS NCs at Hlo « 3.1eV of Schaller et alPand Ellingson 
et alP Error bars show approximate 95% confidence intervals 
reflecting uncertainties due to noise in experimental decays. 



D. Comparison with bulk CM 

In this section we seek to establish a basis for com- 
parison of CM yields between NC samples of different 
sizes and with the bulk material. It has been com- 
mon to compare CM yields at the same scaled ener- 
gies fmj/Eg. This practice follows precedent from the 
bulk impact ionization literature and is useful when con- 
sidering device applications. However, aside from pro- 
viding a convenient way to show data from different 
materials on a single plot, the physical basis for such 
comparisons is not obvious. It may not transparently 
lead to answers of some basic questions, like whether 
or not nanoscale-specific phenomena have a large effect 
on CM. In general, the CM yield for a material sys- 
tem (for example, CdSe or PbS) is expected to be de- 
termined both by particle size and the photon energy, 
y cm (r,Hu)), which can be recast as y crn (Eg,hiu), where 
Eg is the size-dependent bandgap. Much of the existing 
NC CM literature infers an important role for nanoscale 
physics from the fact that their estimates of y cm are 
much larger than reports for y^ k when compared at the 
same relative energy hw/E g . This assumes that with- 
out enhancement y cm (E g ,huj) = y b £ lk [(E bu lk /E g )huj\ , 
or, in other words, that at a given fku, NCs would ex- 
hibit only the CM that would be present in the bulk at 
the lower photon energy (E bulk / E g )hui. To our under- 
standing, though, the only property of y cm (E g ,huj) that 

?bulk , 



a priori scales with E /E g is the energy conservation 

for Huj < 2E g , but this 



requirement, y cm (E g , huj) 
does not seem sufficient to justify the assumption that 
y cm (E g ,huj) = y b c ^ k [(E b g ulk /Eg)fmj] in general as an ad- 
equate description of CM physics in the absence of NC 
enhancement. 
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FIG. 6: (Color online). Diagram of relevant features and 
processes for bulk and NC carrier multiplication. The smooth 
curves are schematics of the bulk lelh and 2e2h densities of 
state, corresponding to X and BX states in an NC. Shown 
for the case of a NC are the lowest X and BX states at E g 
(> _Eg Ulk ) and 2E g and a representative X state formed imme- 
diately after absorption of a high energy photon (hu> 3> 2E g ) 
subject to subsequent intraband relaxation down the X man- 
ifold or coulomb coupling to the BX states. 

To construct a more appropriate reference for compar- 
ison with NC results, we consider bulk material physics 
and explore how y crn (r, hu) behaves if all phenomena ex- 
clusive to the nanoscale are neglected. In the bulk limit 
it is intuitively clear that y cm (r,fkv) is independent of r. 
The competing processes of intraband relaxation and im- 
pact ionization have the same rates for crystals of, say, 
1/im and 0.5 pm, resulting in the same CM efficiency. 
To understand why the impact ionization rate remains 
constant one can start from the first-order perturbation 
theory formulation: 



helh^2e2h = ~T I (2e2/l| V covi | lelh) \ 2 p 2 e2h (E) 

Where V cou \ is the coulomb interaction and p2 e 2h(E) 
is the density of two electron, two hole states (corre- 
sponding to a BX) at the energy E of the initial one 
electron, one hole configuration (which corresponds to 
X in an NC). A reduction in volume has two effects. 
First, the average coulomb coupling is enhanced, with 
|(2e2/i|Koui|lel/i)| 2 oc F" 4 P However, this is fully 
balanced by the reduction in average density of states 
(DOS), since P2e2h(E) oc V 4 . If no new physics are in- 
troduced, this process of shrinking the bulk can be con- 
tinued into the nanoscale with the important conclusion 
that for hoj above the 2E g (r) energy conserving thresh- 
old, y cm (r 7 huj) — y^ k (huj). It must be kept in mind 
that even though the spacing between energy levels is cer- 
tainly larger in NCs, the DOS averaged over sufficiently 
wide intervals is the same as in the bulk, with volume 
scalings px(E) oc V 2 and pbx(E) oc V a . We should 
only expect a divergence from y cm (E g ,huj) = y h c ^ l k {hhj) 



if new physics appear in the nanoscale that have a strong 
influence on the CM process. 

This suggests that comparisons between NC and the 
bulk should be made on an absolute photon energy ba- 
sis as long as fko is well above the energy-conserving 
limit. Then the difference y cm {E gi hui) — y^ k {hJ) would 
be attributable specifically to nanoscale phenomena. In 
contrast, the usual literature comparison at fixed hu/E g 
can significantly exaggerate enhancement over the bulk 
simply because E g (r) > E g ulk , so that, for instance, 
even without novel NC physics, PbSe and PbS NCs with 
E g (r) > 2E g ulk will appear to show at least a two- fold 
CM threshold reduction. It is noted that from a practi- 
cal perspective, bulk- like CM in NCs does indeed present 
a real advantage because the extra carriers can be ex- 
tracted at a higher voltage difference E g (r). Also, a 
hcu/E g basis is useful in comparing a sample's actual CM 
to the maximum possible imposed by energy conserva- 
tion. However, it is not obvious that a comparison of how 
near two different samples are to their separate energy- 
conserving limits can usefully inform on differences in 
their underlying physics. For that, we argue that the 
absolute photon energy basis appears to be more appro- 
priate. 

In light of these considerations, we show in Fig. [7^ 
a summary of literature data on PbS and PbSe NOP*^ 
rcplotted on an absolute energy axis along with values 
of CM yields in bulk PbS films reported by Smith and 
DuttonpS These authors studied the photoconductivity 
of commercial PbS films and found an increase in pho- 
tocurrent response at shorter wavelengths which they at- 
tributed to a CM process, emerging from a threshold 
huo « 2eV and rising approximately linearly to j/J™ 1 w 2 
at hu> — 6eV. It should be kept in mind that there are 
numerous potential sources of error in this bulk CM de- 
termination, some of which we detail later, but it is nev- 
ertheless interesting to note that the CM yields for NCs 
reported in the literature appear only modestly enhanced 
over these bulk values. Except in the case of very large 
(small E g ) NCs following the universal curve of Schaller 
et al., CM yields are within a factor of « 2 of the bulk re- 
port, and exhibit a similar CM energy threshold between 
2-3eV. Figure [TJd shows our own estimates compared to 
the NC and bulk literature data at hu> — 2>.leV plotted 
against E g (r). Our results are below the bulk CM re- 
ported value, those of Ellingson ct aP appear consistent 
with it, and the Schaller et al. 5 results fall well above for 
larger NCs. 

Reaching a robust conclusion at this stage on the rel- 
ative strengths of CM in bulk and NC forms is difficult 
because of potential uncertainties in the bulk values re- 
ported Smith and Dutton. First, the authors did not 
present a characterization of the commercial PbS films 
studied, and it is possible that significant oxidation may 
have taken place since no protective coating was usedJ^ 
This is important since exposure to O2 is known to cause 
significant changes in bulk PbS photoconductivity.^ Sec- 
ond, the reported yields are very sensitive to any sys- 
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FIG. 7: (Color online) (a) Summary of CM reports in the 
literature for PbSe and PbS NCs by Schaller et alPand Elling- 
son et alP" compared on an absolute excitation e nerg y scale 
to bulk PbS values reported by Smith and Dutton! — ' For the 
Ellingson et al. data, filled circles or open shapes indicate 
data taken with intraband or band-edge probe beams, while 
the colors blue, green, red, and black denote E g = 0.72cV, 
0.82eV, 0.9feV PbSe and 0.85eV PbS NC samples. 11 Two 
curves are displayed based on the Schaller et al. universal 
huj/E g dependence for E g = 0.6eV and E g = 0.9eV. (b) 
Plot of CM yields against the size dependent bandgap E g for 
huj — 3.1eV. The dashed line is the CM yield reported for the 
bulk at this excitation energy. 



tematic errors in determining the number of photons ab- 
sorbed by the film. The third complication is the possible 
variation of photoconductive gain with Hlo. For exam- 
ple, at blue wavelengths carriers are generated on average 
closer to the film surface, where the greatest concentra- 
tion of trap states are expected to reside. Moreover, it is 
difficult to say a priori whether the gain would increase 
or decrease. These considerations highlight the need for 
a careful determination of CM in bulk films of PbS and 
PbSe before a definitive comparison with NCs can be 
made. With the data at present it is difficult to conclude 



that nanoscale phenomena are responsible for strong CM 
enhancement, as we have discussed in the previous para- 
graph and in Fig. [7] 

Given the possibility that CM in NCs might follow 
largely bulk- like physics, it is interesting to examine what 
is known about the NC-specific physical mechanisms that 
could affect the multiplication process. The most com- 
monly cited rationalizations of CM enhancement in NCs 
are the possibility of strong coulomb interaction and slow 
intraband relaxation.^ It could be argued, for exam- 
ple, that coulomb couplings are significantly enhanced in 
the nanoscale based on the much faster Auger relaxation 
rates of band-edge multiexcitons compared to the bulk. 
This enhancement of Auger rates at the band edge is 
thought to be due to a relaxation of momentum conserva- 
tion requirements brou ght a bout by the finite nature and 
abrupt surface of NCsJS2E21 However, because momentum 
conservation is not a limiting constraint on impact ioniza- 
tion in the bulk at high excess kinetic energies,^ and it 
is not clear that the nanocrystalline form should still ex- 
hibit significant enhancement. Calculations by Allan and 
Delerue suggest that fci e i/i— >2e2h is if anything smaller in 
PbSe NCs than for the bulkP 

Similarly, there is still no evidence of a phonon bottle- 
neck for intraband relaxation at high electron and hole 
kinetic energies. Due to practical considerations relating 
to experimental time resolution, most studies on NCs 
have focused only on relaxation from so me of the lowest 
excited states to the band edgeP^iM] Even then, they 
find very fast picosecond relaxation times. Moreover, at 
the high excess kinetic energies required for CM, the X 
manifold is much denser and it seems less likely that a 
phonon bottleneck effect could play a very large role. 

The remaining potential nanoscale CM enhancement 
mechanisms have to do with the discrete state structure. 
Certainly, the discrete nature of states in a NC is critical 
near the energy conservation threshold, as no CM can 
occur when hu; < 2E g (r) even though the bulk 2e-2h 
DOS is finite. However, if we restrict our attention to 
hco well over 2E g (r), as has been the case when large 
CM yields have been reportedP^ it is plausible that the 
BX manifold is sufficiently dense that bulk-like behavior 
could result. Further, even if there were deviations, we 
would not expect them to be monotonic in either E g (r) 
or hut. Finally, it is possible that there could be strong 
coupling between X and BX,ES3 but not enough is known 
about phase and population relaxation mechanisms of 
carriers with high kinetic energies to conclude that such 
effects would be important for CM. 

All these arguments above should not be taken as proof 
or justification that y cm (E g , huj) = y^ k (hcu) for NCs, 
but simply to show that such a conclusion would not be 
inconsistent with what is experimentally known about 
NC photophysics. Too little is understood about the 
physics of highly energetic carriers in NCs to make strong 
a priori predictions of the role of nanoscale phenomena 
in CM. 
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IV. CONCLUSIONS 

The principal experimental conclusion of this work is 
that CM yields in our PbSe and PbS NC samples esti- 
mated by transient photolumincscence are well below the 
values that have been reported in the literature for PbSe 
and PbS NCs using transient pump-probe techniques!^ 
It should be noted that these previous reports them- 
selves show significant numerical disagreement between 
each other even though they employ nominally similar 
methods. In broad terms, the variation between the re- 
ports of Schaller et al., Ellingson et al., and our own must 
ultimately stem from either systematic differences in data 
acquisition procedures, variation in the way CM is deter- 
mined from observed decays, or actual sample-to-samplc 
differences of the CM efficiency. The fact that Ellingson 
et al. and Schaller et al. use nearly equivalent methods 
for estimating CM but find conflicting results suggests 
that their samples are inherently different, or that the 
two groups handled these samples differently during the 
course of their experiments. In our own work there is a 
possibility of systematic error related to the calibration 
method we use, but we have argued above that this alone 
cannot readily account for the contrast with the existing 
literature. The answer may yet lie in sample-to-sample 
CM variation, and if so would suggest that CM in NCs 
is strongly affected by defects or surface ligand type and 
coverage. 

The second effort of this work has been to establish a 
basis for comparing CM efficiencies in NCs and the bulk 
that more clearly isolates the effects of changes in under- 
lying physics. We have argued that an absolute photon 
energy basis is more appropriate than the usual hw/E g 
approach for this purpose, and by comparison to values 
reported for the bulk, we found that the CM yields re- 
ported for NCs do not immediately suggest a very large 
role for nanoscale-specific phenomena. Because these 
bulk values themselves could be beset by large errors, 
it is difficult to reach a definite conclusion. A modern, 
robust, assessment of CM in bulk PbS and PbSe will be 
necessary for this to be possible. Similarly, understand- 
ing the variation in the NC CM literature will require 
applying multiple experimental methodologies to identi- 
cal NC samples or, more importantly, the development of 
new spectroscopic techniques that are more specifically 
taylored to multiexciton detection than the population- 
dynamics based measurements in use at this time. A 
clear picture of the CM process in the transition from 
the bulk to the nanoscale will have to wait for experi- 
mental efforts on these two fronts. 
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APPENDIX A: POPULATION MODELING 

As has been described in detail in our previous work on 
CdSe NCs the exponential components in a tPL decay 
can be related to the MX and X populations immediately 
following excitation through the following approximate 
expressions: 

a x oc k r x d V\ 
a BX oc (k r B a x - k r x ad ) P2 

where p\ and pi are the population of NCs that start 
with at least an exciton or a biexciton respectively at 
time 0. These populations are given by p m = ^2u> m Iki 
where the population of the fc-th multiexciton state, Ik, 
is determined by poisson statistics, taking into account 
excitation beam inhomogeneity and position dependent 
collection efficiency: 

h = J 0(z)^e-"^d 3 x n(x) = aj(x) 

where j{x) and </>(x) are the photon flux and collection 
efficiency at position x, and a is the absorption cross sec- 
tion. n(x) is the average number of photons absorbed per 
pulse by an NC located at x. If j(x), (f>(x) and a were 
known, it would be possible to compute the Ik up to a 
common proportionality constant and obtain, by compar- 
ison with experiment, the saturation ratio (asx /a>x)sati 
the value of asx/o-x when p 2 — Pi- However, both j(x) 
and especially 4>{x) are difficult to determine accurately 
in our apparatus. Instead, we exploit the fact that the 
shape of p\ as a function of excitation power fully de- 
termines the shape of p 2 - To see this, we note that 
during any of our experimental power series n(x) only 
changes in magnitude while retaining its shape. Setting 
n{x) — noh(x), where h(x) is a fixed shape and no is a 
constant parameter, one can show that: 

P2{n j =Pi{n ) ~n - — 
on 

In the above, no can be replaced with any quantity 
proportional to it, such as average excitation power, so 
knowledge of the absorption cross section is not required. 
Therefore, if one finds any h(x) and <j>(x) so that the 
calculated pi(no) closely fits the shape of the observed 
a x excitation series, then the P2(no) calculated with the 
same h(x) and <fi(x) will be proportional to asx- The 
results of this procedure, shown in Fig. [3J3, demonstrate 
that the a bx we observe match very well the trend we in- 
dependently predict from the ax evolution, further sup- 
porting our assignment of this fast component in the tPL 
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to the biexciton. Our estimate of (asx /&x)sat is then 
obtained as the ratio of the proportionality constants 
relating ax to pi and asx to P2- We find saturation 
values (cibx /ax)sat of 2.5-4 using this method. Since 
{o-bx I 'ax) sat = k r gx/k r x d — 1, the corresponding values 
of k^x /k r x d ar e in the range 3.5-5. 

APPENDIX B: BX AND X RADIATIVE RATES 

We present here a calculation of k r %x I k r x d for a sim- 
ple model of the lead chalcogenide ground state. The 
lS e and ISh states in lead chalcogenide are eight-fold de- 
generate. There are four equivalent valleys in the band 
structure and two-fold spin degeneracy. The possible X 
electronic configurations can be labeled i e m^, and the 
BX configurations i e j e mhnh, where i,j,m,n G 1...8. 
Because total momentum and spin must be conserved 
during an optical interaction, only the recombination of 
an electron and hole with the same k and same spin is 



allowed. Assuming that particle momentum and spin 
remain good quantum numbers, each electron state is 
connected by a dipole transition to exactly one of the 
eight hole states. By symmetry, these transition dipole 
moments all have the same magnitude \(x\. We can then 
calculate the radiative square transition dipole of each X 
and BX microstate. In the case of X, there are 8 con- 
figurations of type l e lh with k ra d — A* 2 ; & n d 8 • 7 of 
type l e 2; l with k ra d = 0. Similarly, for the BX, there 
are ( 2 ) configurations like l e 2 e lh2h with k ra d — fJ- 2 + fJ 2 , 
8-7-6 configurations of type l e 2 e lh2>h with k ra d = A' 2 ; 
and (2) (2) dark l e 2 e 3/ l 4? l -type states. Taking the ther- 
mal average, one finds fc™ = [J, 2 / 8 and k r ^x — A^/2, 
and therefore k r ^x — 4/c™ d . This result should remain 
approximately valid even in the presence of perturbations 
that mix states with different quantum numbers or cou- 
ple the electrons and holes, as long as the width of the 
resulting energy fine structure is sufficiently smaller than 
the available thermal energy, kT. 
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In the literature CM yields are often reported as an internal 
quantum efficiency (IQE), which is related to our y cm by 
IQE = 100%(3/cm + 1). 

As in our previous workp^we delay fitting of the 1.55eV- 
excited decays by a time ~ tbx/2 to minimize unwanted 
interference from higher MX tPL components in the deter- 
minations of ax and asx- 

For example, we find tbx ~ 58ps and tbx ~ 140ps for 
E g = 0.84eV and E g = 0.68eV PbSe NCs respectively, 
while Beard et alP have determined tbx = 67ps for E g = 
0.84eV and Schaller et alf^ report tbx = 149ps for E g = 

0.64eV. 

The steep volume dependence |(2e2/i| V C oui|lel/i}| 2 oc V~ 4 
might appear surprising at first. It should be kept in 
mind that this square matrix element is averaged over 
all 2e2h configurations of nearby energy. However, for 
{2e2h\ Vcouil lelh) ^ 0, conservation of momentum and 
spin must be satisfied and one of either the initial elec- 
tron or hole must not change state. The proportion of fi- 
nal 2e2h states that violate these conditions and thus have 
{2e2h\ y c0 ui| lelh) — increases with volume, ultimately re- 
sulting in a stronger volume scaling of |(2e2/i|Vcoui|lel/i}| 2 
than would be expected from averaging only the non-zero 
terms. 



